Abstract. We introduce the mixing resulting from the 
radiative diffusivity associated with the radiative viscos- 
ity in the calculation of stellar evolution models. We find 
that the radiative diffusivity significantly diminishes the 
efficiency of the gravitational settling in the external lay- 
ers of stellar models corresponding to types earlier than 
w G. The surface abundances of chemical species predicted 
by the models are successfully compared with the abun- 
dances determined in members of the Hyades open cluster. 
£N| ' Our modeling depends on an efficiency parameter, which 
. is evaluated to a value close to unity, that we calibrate in 
this study. 
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1. Introduction 

Microscopic diffusion, sometimes named "atomic" or "ele- 
ment" diffusion, when used in the computation of models 



of main sequence stars produces (see Chaboyer et al. 2001 
for an abridged revue) a change in the surface abundances 
from their primordial values and a shift of the evolution- 
ary tracks towards lower luminosities and temperatures 
(Morel & Baglin 1999| ), consequently increasing the ex- 
pected age of binaries and clusters. For the Sun, the inser- 
tion of microscopic diffusion has spectacularly improved 
the agreement between the theoretical and t he "seismic" 
model (e.g. Christensen-Dalsgaard et al. 1996 ). Indeed, for 
low mass stars and for population II stars, atomic diffusion 
eliminates, at least partially, discrepancies betw een ob ser- 
vations and theor etical models (Mor el & Baglin 1999; Le- 



breton et al. |1999[ Ca yrel et al. |1999| ; Chaboyer et al. |2001 
Salaris & Weiss |200l| ). 

For stellar models with mass larger than w 1.4 Mq, the 
use of microscopic diffusion alone produces at the surface 
an important depletion of helium and heavy elements and 
a concomitant enhancement of the hydrogen content. As 
an example, for a 2 Mq model computed with atomic dif- 
fusion, at the age of 20 Myr the surface mass fraction of hy- 
drogen increased from X = 0.70 to X (Z 0.98, while helium 
and heavy elements are strongly depleted. This strong de- 
pletion, which increases with the mass of th e star, is not 
observed either by Varenne & Monier (1999) in stars be- 
longing to the Hyades open cluster or in OB associations 
(Daflon et al. 2000) and in the Orion association (Cunha 
& Lambert [1994D . 

There is no possibility of observing helium in main se- 
quence stars of intermediate mass to check its depletion; 
such stars do not show low surface metallicities. This large 
helium depletion is not likely to be observed in main se- 
quence B-stars, as supported by their spectral classifica- 
tion which is mainly based on the relative strength of the 
helium spectral lines. 

The large efficiency of the atomic diffusion in the outer 
layers results from the large temperature and pressure gra- 



dients. For models of stars with types later than ss G, 
i.e. M* 5s> 1,2 Mq, the external convection zone acts as 
a "reservoir" which inhibits atomic diffusion. On the con- 
trary, for stars with masses larger than <1 1.4 Mq, the ex- 
ternal helium and hydrogen convection zones are so tiny 
that the reservoir is not large enough to be efficient. 

Such misleading physics in the calculation of micro- 



scopic diffusion coefficients (e.g. Cox et al. 1989; Proffitt & 
Michaud|l~99ll; Michaud & Proffit |1993|; Thoul et al. |1994[) 



are very unlikely owing to the high degree of sophistica- 
tion of the method based on the kinematical theory of 



gases (Burgers 1969) and also as seen by its success when 
applied to the Sun. Clearly, a physical process which in- 
hibits the efficiency of microscopic diffusion is lacking. 

Among physical processes acting against gravity, re- 
cent developments were focused on: 

— The radiative accelerations, which are efficient only in 
the case of chemical species with a large charge, e.g. 
iron. As such "metals" have small relative abundances, 
the radiative accelerations are not efficient limiting the 
sedimentation of helium and heavy elements as a whole 
(e.g. Alecian et al. |l989j ; Turcotte et al. |l998a ). 

— The turbulent mixing generated by the shear created 
by the differential rotation between the solar convec- 



tion zone and the radiative core (Richard et al. 1996 



Gabriel 1997; Brun et al. 1998). It has been recently 
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claimed (Schatzman et al. 200C) that the shear insta- 
bility invoked to produce this turbulent mixing only 
concerns a narrow region and may be not as efficient 
as expected. 

The rotation-induced mixing accounting for the trans- 
port of angular momentum explains well the hot side of 
the Li dip in the Hyades (Talon & Charbonncl [1998 ). 
The use of a hypothetical mixed reservoir with an ad- 
hoc adjusted size (e.g. Turcotte et al. [1998b] ). 
The refinement of the physical description. In the cal- 
culations of G to A type star models (Turcotte et 
al. [1998a] , |1998b| ; Richer et al. |2000| ), the local mean 
Rosseland opacity is updated at each time step, ac- 
cording to local changes in the density and tempera- 
ture, but also in the genuine mass fractions of chemical 
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species. In some favorable situations the authors have 
obtained the formation of an "iron" mixed convective 
zone that inhibits the efficiency of the atomic diffusion. 
Even with this noticeable improvement of the model- 
ing, several difficulties remain; in particular the mod- 
els do not reproduce the abundances of star members 
of the Hyades open cluster as observed by Varenne & 



Monier (1999); their figures 5 to 8 for elements from 
C to Ba only reveal a mild depletion by a factor w 3 
for 6 600 K< T cff < 7 800 K. Also the the oretica l strong- 
depletions predicted by Turcotte et al. ( 1998b ) for ef- 
fective temperatures larger than 6 600 K are not ob- 
served. 

— The influence of the stellar mass -loss described phe- 
nomenologically (Chaboyer et al. 1999] ) in the model- 
ing of Procyon A. 

In this exploratory work, we focus on the role of the 
radiative diffusivity generated by the photon-ion collisions 
that is not presently taken into account in the microscopic 
diffusion coefficients. As we shall see, the photon-ion col- 
lisions are an efficient physical process that inhibits the 
large sedimentation of helium and heavy elements in outer 
layers of main sequence star models with types earlier than 

In Sect. H we present a simplified phenomenological 
model of the radiative diffusivity resulting from photon- 
ion collisions. Section || briefly describes the theoretical 
framework used for constructing stellar evolution models. 
In Sect. ^ we estimate the radiative diffusivity parame- 
ter connecting the radiative kinetic viscosity to the ra- 
diative diffusivity; our investigation is based mainly on 
observed abundances of C, O, Fe and Mg in stars belong- 
ing to the Hyades open cluster as reported by Varenne & 
Monier (1999). Illustrations are given in Sect. || and, in 
Sect. ^, we summarize and discuss the main results of the 
present investigation. 



2. Radiative diffusivity 

Apparently, the radiative diffusivity, first used by 



Baglin (1972), has not been introduced for intermediate 
stellar mass models, although implemented in the evolu- 
tion of models of larger mass (Heger et al. 2000). Accord- 
ing to the standard description of the micro-physics of 



gases (e.g. Burgers |1969| ; Mihalas & Weibel-Mihalas |1984 , 
Chap. 1) the viscous stresses due to electric forces tie up 
the constituent particles, electrons and ions. In an astro- 
physical plasma with radiation, the kinematic viscosity v 
has two components: 

V = VjnoX + ^rad • 

The kinematic molecular viscosity, t'moi, results from en- 
ergy exchanges between thermal collisions leading to exci- 
tation and ionization of atoms and ions. Thus, the thermal 
collisions generate a "molecular mixing" resulting from the 
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Fig. 1. Profiles of kinetic radiative viscosity i/ ra( j (heavy 
full), kinetic molecular viscosity v mo \ (heavy dashed), hy- 
drogen diffusivity c?h (dot-dash-dot-dash) and helium dif- 
fusivity due (dotted) with respect to radius in the radia- 
tive zones of a IMq (R < 0.61 Rq) and in a 1.8 Mq 
stellar models at age 700 Myr. The kinematic molecular 
and radiative viscosities and the diffusivities are in c.g.s. 
(cm 2 s _1 ) units. 



molecular diffusivity associated with v mo \. The kinematic 



molecular viscosity is approximated by (Schatzman 1977): 

^ mol ^2.21Hr 15 — i±Z£, (i) 
p 3mA 

for a mixture of hydrogen and helium, assuming that the 
heavy elements have a negligible effect. T is the tempera- 
ture, p is the density, X is the hydrogen mass fraction and 
In A is the so-called Coulomb logarithm (e.g. Michaud & 
Proffit |1993| ). With the more elaborate description based 
on the kinetic theory of gases (Burgers 1969), the molecu- 
lar diffusion has two components: the isotropic one, being 
the molecular diffusivity that generates a mixing and the 
anisotropic one, responsible for the gravitational sedimen- 
tation in stars. 

The kinematic radiative viscosity j/ ra( j arises because 
photons deposit their momentum in the fluid element into 
which they are absorbed. Thus, the radiative collisions 
generate a "radiative mixing" resulting from the radia- 
tive diffusivity associated with t'rad- In the optically thick 
medium the kinematic radiative viscosity is expressed as 
(Thomas |1930| ; Mihalas & Weibel-Mihalas [1984 P- 461- 
472): 

4 «T 4 

^=15^' (2) 
where a is the radiation density constant, c is the speed 
of light in vacuum and k is the mean Rosseland opacity. 
As already emphasized by several authors more than 



seventy years ago (e.g. Milne 1930), the radiative viscos- 



ity is large in the external layers of stars (Baglin 1972 ). 
Figure [l] shows the profiles, with respect to radius, of v mo \ 
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and I'rad in the radiative zones of two main sequence stellar 
models of respectively 1.0 M Q and 1.8 M at age 700 My. 
In the radiative zone of the 1M Q model (R < 0.61 Rq), 
the total kinematic viscosity mainly results from the kine- 
matic molecular viscosity while in the radiative outermost 
layers of the 1.8 Mq model, the kinematic radiative viscos- 
ity overcomes the kinematic molecular viscosity by sev- 
eral orders of magnitude. The large value of the radia- 
tive viscosity in outer layers of stars is associated with a 
large radiative diffusivity that makes a significant contri- 
bution to the mixing. Like the diffusion due to molecular 
collisions, the radiative diffusion has an isotropic compo- 
nent, the radiative diffusivity generating mixing, and an 
anisotropic component causing the radiative accelerations 
that act against gravity (Alccian 1996 ). 

For stars with mass Ss 1.2M© the effect of radia- 
tive diffusion is swallowed up by the mixing in the ex- 
ternal convection zone, as in the Sun. For stars of types 
earlier than G, with increasing mass (i.e. effective temper- 
ature) the external convective zones become smaller and, 
in outer layers, the radiative diffusivity is a process which 
acts more and more efficiently against the gravitational 
sedimentation. 

In the most elaborated stellar evolution codes, the co- 
efficients of microscopic diffusion are derived using the 
Burgers (|1969[ ) formalism (e.g. Coxet al. |1989| ; Pro ffitt fc 
Michaud 1991; Michaud & Proffit [199$ Thoul et al. |1994fl . 
But, in the standard framework of the calculation of col- 



lision integrals and r esistance coefficients (Burgers 1969 ; 
Paquette et al. 1986 ), only the ion-ion and electron-ion 
collisions are taken into account and the derived diffusiv- 
ity corresponds only to v mo \. 

Figure |l| also shows the profiles of the hydrogen c?h 
and helium due molecular diffusivities. Note that (in, (fee 
and v mo \ have similar profiles. This behavior suggests a 
linear relationship, with a dimensionless coefficient close 
to unity between the molecular diffusivity and the kinetic 
molecular viscosity. 

It is out of the scope of this exploratory work to com- 
pute the diffusion coefficients including both the ther- 
mal a nd the radiative collisions according to the Burgers' 
(|l969j Sect. 18 & 60) formalis m. We express the diffusiv- 
ity associated whith z/ ra( j by a diffusion coefficient propor- 
tional to the kinematic radiative viscosity: 

rfrad = Dr t'rad- 

The efficiency factor, Dr, is a dimensionless phenomeno- 
logical parameter to be calibrated. Hereafter, we shall refer 
to Dr as the "radiative diffusivity parameter" . We assume 
that Dr is constant even if, a priori, there is no reason 
for a unique value to apply for all atoms and ions of all 
chemical species and in all stellar conditions. 

To illustrate the effect of the radiative diffusivity, using 
the results of Sect. ||, Fig. || shows the abundance profiles 
with respect to radius in a 1.8 Mq main sequence stellar 
model at age 570 Myr computed respectively with (Dr = 



1.0) and without (Dr = 0.0) radiative diffusivity. Note the 
small surface abundances of metals predicted by the model 
computed without radiative diffusivity and the realistic 
values obtained with Dr = 1.0. 

The simplified description of the radiative diffusivity 
introduced in this paper is phenomenological and prelim- 
inary. Particularly it ignores the ionization stage of ele- 
ments which could modify the efficiency of the radiative 
diffusivity of chemical species under consideration. It ob- 
viously needs further extended development based on the 
kinetic theory of gases with radiation. 

3. Computation of models and physical inputs 

Stellar modeling: Stellar models have been computed us- 
ing the CESAM code (Morel 1997). Th e physics employed 
is the same as in Morel et al. ( 2000 ). The stellar evolu- 
tion sequences are initialized with homogeneous pre-main 
sequence models. We follow in detail the abundances of 
all chemical species of our diffusion network. The set of 
evolutionary differential equations fulfilled by the abun- 
dances of the isotopes is numerically integrated using a 
finite element method^]. The microscopic diffusion coeffi- 
cients are derived from the set of Burgers' ( 1969] , formula 
18.1-18.2) flow equations. The charge of a given isotope is 
taken as the average^ charge over all its ionization states. 
The collision integrals and resistance coefficients are taken 
from the tabulations of Paquette et al. ( |1986| ). For each 
ion, the coefficient of radiative diffusivity is added to the 
coefficient of molecular diffusivity. 

Input physics: The eff equation of state (Eggleton et al. 



1973) was used, together with opacities determined by 
Rogers et al. (|1996| ) using the solar mixture of Grevesse 



& Noels (1993), complemented at low temperatures by 
Alexander & Ferguson (1994) opacities. In the convection 
zones the mixing is insured by a large turbulent diffusion 
coefficient and the temperature grad ient is computed ac- 
cording to the Bohm-Vitense ( 1958 ) convection theory; 
the mixing-length parameter is fixed to the value in the 
standard solar model computed with the same physics^ 
The atm osphe re is restored using Hopf 's atmospheric law 
(Mihalas 1978 ). The element diffusion is not computed in 
the atmosphere where the chemical composition is fixed 
at the value in the outermost shell of the envelope. 

Nuclear network: The relevant nuclear reaction rates are 



taken from the NACRE compilation (Angulo et al. 1999 ). 
The nuclear network contains the isotopes: 1 H, 2 H, 3 He, 
4 He, 7 Li, 7 Be, 12 C, 13 C, 14 N, 15 N, 16 and 17 0. They are 



1 The derivation of the equations and the algorithms are de- 
tailed in the "Notice de CESAM" and available on request o n 



the WEB site http://www.obs-nice.fr/morel/CESAM.html 

2 Weighted by the ionization rates. 

3 But the inefficient radiative diffusivity, see Sect. W. 
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Fig. 2. Profiles, with respect to radius, of abundance mass fractions in a 1.8 Mq "Hyades" stellar model at age 570 Myr 
computed with Dr = (left) and Z?r = 1 (right). H x 1 (heavy full), He x 1 (dashed), C x 100 (dot-dash-dot-dash), 
N x 100 (dotted), O x 100 (dash-dot-dot-dot-dash), Fe x 100 (thin full). 



processed by the PP+CNO thermonuclear reactions rele- 
vant for the main sequence evolution. The elements 2 H 
and 7 Be are set at equilibrium. 

Diffusion network: Along the main sequence the chemi- 
cal species, heavier than O, are not nuclearly processed. 
However, their abundances vary with respect to time ow- 
ing to diffusion processes, contributing to changes of Z, 
thus generating opacity variations. Morel et al. (1997) 
introduced a dummy "mean metal" to approximate the 
mean opacity changes caused by the variations of Z re- 
sulting from microscopic diffusion. To mimic more closely 
the effects of low and high charged chemical species, we 
use here two dummy mean metals. The first simulates the 
metals heavier than potassium. For the solar mixture of 
Grevesse & Noels (1993) used so far, their mearj0 charge 
and mean atomic weight are respectively 26 and 55. There- 
fore, a dummy "iron" , labeled hereafter { 55 Fe}, is well rep- 
resentative of the heaviest chemical species. Indeed, the 
second dummy metal simulates the isotopes with atomic 
numbers from 9 to 20, and also complements the mix- 
ture. Its charge and atomic mass are respectively 12 and 
24, corresponding to a "dummy" magnesium { 24 Mg}. For 
the sake of simplicity we shall approximate { 55 Fe} and 
{ 24 Mg} with respectively the iron and the magnesium 
abundances measured spectroscopically. 



Initial abundances: To fulfill the basic relationship X{ n i + 
Yini + Zj n i = 1, the initial abundances of isotopes are de- 
rived (1) from the isotopic ratios of nuclides (Anders & 
Grevesse |l989|) , (2) fro m the initial mass fraction of he- 
lium, Yi n i = 3 Hcj n i + 4 Hci n i, and (3) from the initial mass 
fraction of heavy elements relative to hydrogen (-j-)mi- 
As deuteron is assumed at equilibrium the initial amount 



of 3 He is enhanced by the initial meteoritic amount of 
deuteron (Geiss & Reeves 1972 ). 

Failing anything better, we use the meteoritic val- 
ues of Anders & Grevesse ( 1989 ) for the initial fractions, 
by number, of metals within Z: 12 C + 13 C = 0.2455, 
i4 N+ is N = 0.06458, 16 0+ 17 = 0.5130. From these data 
one infers the initial abundances for the dummy metals 
{ 24 Mg} = 0.1631, { 55 Fe} = 0.01563. For the initial iso- 
topic ratio (by number) we use 3 He/ 4 He= 1.1 10~ 4 (Gau- 
tier & Morel |1997| ). 



Weighted by number. 



For any chemical species X we shall refer to [^] as the 
quantity: 

[f)^(f) t - lo6 (£L' < 3 » 

which is a measurement of the depletion of X at the age 
t, under consideration. 



4. Constrains on the radiative diffusivity 
parameter. 

The radiative diffusivity parameter Z?r is estimated by 
comparing relevant observations to theoretical results. 
Firstly, owing to the high accuracy of helioseismological 
observations, some insights can be gained using the stan- 
dard solar model despite the inhibition, by the convective 
mixing, of the microscopic diffusion in the outer layers. 
Secondly, the most sensitive observations, relevant to the 
atomic diffusion, are the surface abundances of chemical 
species of members of well-observed open clusters. When 
the age of the cluster is not older than « 1 Gyr, the effects 
of element diffusion remain small for k, G dwarfs and the 
initial abundances can be reasonably inferred. For inter- 
mediate mass stars we shall estimate the value of the ra- 
diative diffusivity parameter using the observational ma- 
terial of the Hyades open cluster and, for early type stars, 
we shall use observations of Orion association. 
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Table 1. Depth of the convection zone (solar unit), sur- 
face isotopic ratio of 3 He versus 4 He, surface helium mass 
fraction Y s , small difference Svq2 (a*Hz) and [Li] s lithium 
abundance in dex (H = 12) in calibrated solar models 
computed with different values for the radiative diffusiv- 
ity parameter Dr. 



D R 


-Rzc 






8VQ2 


[Li] s 





0.7127 


4.343 


0.2430 


9.09 


2.603 


1 


0.7128 


4.345 


0.2430 


9.10 


2.609 


5 


0.7129 


4.353 


0.2431 


9.12 


2.608 


10 


0.7131 


4.359 


0.2433 


9.15 


2.607 


50 


0.7142 


4.381 


0.2448 


9.32 


2.502 



Observations 
J? zc = 0.713 ± 0.001 
[^g£] a = 4.4 ± 0.4 
Y s = 0.240 ± 0.009 
Svw = 9.01 ±0.10 



References 
Basu & Antia (|l995| ) 
Bodmer et al. (1995) 



Basu ( |l997j ) 

Grec et al. ( |t997|) 



Solar models: In the external layers of a solar model, the 
radiative diffusivity is swallowed up by the convective mix- 
ing. Beneath the convection zone, where precise informa- 
tions on the internal structure result from helioseismology, 
the effects of radiative diffusivity should be compatible 
with the relevant properties of the nowadays very precise 
standard solar model. This assertion leads to an estimate 
of an upper limit for Dr. 

We have computed calibrated solar models for the fol- 
lowing values of the radiative diffusion parameter Dr = 
0, 1, 5, 10, 50. To compare solar models with helioseismo- 
logical observations, we use improved physical data and 
standard solar initial chemical composition (details can 
be found in Morel et al. 1997). In the range 0.1 i? ^ R ^ 
0.9 Rq where the inversions of the helioseismic data are 
reliable, Fig. || shows the relative differences between the 
sound speed in the calibrated solar models and the seis- 
mic s ound speed "experimental" results of Turck-Chieze 
et al. ( 19971 ). F° r the models computed with Dr ^ 10, the 
relative discrepancy between the sound speed in the Sun 
and in the model are not larger than a few 10~ 3 . Table |] 
shows that the theoretical values of standard solar charac- 
teristics are all compatible with observations; nonetheless, 
these data are not strongly dependant on Dr. 

From Fig. |^, we estimate that a reasonable upper limit 
for the radiative diffusivity parameter is close to [Dr] ^ 
10. 
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Fig. 3. Relative differences in sound velocities between 
the Sun and calibrated solar models with Dr = 
0, 1, 5, 10, 50 respectively: full, dashed, dot-dash-dot- 
dash, dotted, dash-dot-dot-dot-dash. 



spect to the effective temperature, i.e. masses, neither for 
C, N and O nor for heavier elements, as it is shown in 
Fig. 3 of D aflon et al. fl2000| ) or in Fig. 17 of Cunha & Lam- 
bert (1994). These recent star- forming regions are younger 
than 20Myr and most of the stars analysed in Cunha & 
Lambert (1994) have an age supposed to be around 7Myr. 

In models of such massive stars including element diffu- 
sion alone (Dr = 0.), helium and metals empty the outer 
layers in an evolutionary time shorter than ss 10 Myr. In 
these models the outer hydrogen and helium convection 
zones are so tiny that the convective mixing is ineffec- 
tive. Resulting from the large pressure and temperature 
gradients, the element diffusion is very efficient, leading 
to a large sedimentation of helium and heavy elements. 
The inclusion of the radiative diffusivity generates a mix- 
ing which avoids this unrealistic result. Due to the large 
effective temperature (T e ff 2b 10 000 K), the kinematic ra- 
diative viscosity is large in the outer layers, see Eq. (||), 
and so the radiative diffusivity; that generates an efficient 
mixing inhibiting the large efficiency of the element diffu- 
sion and allows one to reproduce the observed absence of 
depletion of helium and metals. 

To estimate Dm, adopting a mean age of 7 Myr, we 
have computed the evolved abundances of C, N, O and 
Fe for stellar masses from 7 to 17M Q for different values 
of the radiative diffusivity parameter. In these early type 
star models the radiative diffusivity inhibits the efficiency 
of the element diffusion as soon as Dr is larger than the 
lower limit |D R J = 0.05. 



Massive stars: Massive stars, belonging to OB associa- 
tions, provide a test of efficiency for the radiative diffusiv- 
ity. For stars having masses from 7 to 17 Mq observations 
do not reveal any chemical abundances depletion with re- 



Hyades: The Hyades cluster is a nice laboratory for test- 
ing stellar structure and evolution theory because of its 
numerous well analysed stars and its well measured dis- 
tance. We have computed star models taking into account 
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the atomic diffusion including radiative diffusivity for 
different values of Dr. We refer to Lebreton et al. (2001) 
for the cluster parameters. We adopt a metallicity of 
= +0.14 ± 0.03 dex and an age of 550 ± 50 Myr as our 
models are computed without overshoot of the convective 
core. Alike for the Sun, a typical G dwarf Hyades member 
of 1.1 Mq is not affected by the radiative diffusivity but, 
owing to the element diffusion, the initial helium mass 
fraction needs to be enhanced up to Yj n i = 0.259 to 
ends at Y — 0.255 after 550 Myr of evolution. Also, the 
initial metallicity needs to be [Fe/H]i„i = 0.147 dex to 
reach [Fe/H] = 0.140 dex at the age of Hyades. For our 
"Hyades" models we use the following initial mass frac- 
tions (in dex with respect to H e 12): 3 He=8.056, 
4 He=11.56, 12 C=9.738, 13 C=7.818, 14 N=9.299, 
15 N=6.896, 16 O = 10.24, 17 0=6.847, { 24 Mg}=9.416, 
{ 55 Fe}=9.838. 

Figure || shows the surface abundances of C, O, Mg, 
and Fe with respect to effective temperature in mod- 
els computed with Dr = 0.0, 0.5, 1.0, 5.0, 9. 0, tog ether 
with the observations of Varenne & Monier (1999). For 
main sequence stars, with effective temperatures between 
6 600K< T c ff < 7 800K, the depletions of abundances are 
decreasing with respect to the radiative diffusivity param- 
eter. With Dr = 1 the depletions amount to about a factor 
« 2 — 3 for all chemical species; that corresponds to the 
observations even if, for oxygen, the distribution of the ob- 
served abundances are rather scattered, probably resulting 
from non-LTE effects or 3D temperature inhomogeneities 



in the atmosphere (Asplund et al. 2000a). 

Taking into account observational error bars, Fig. |] 
shows that the best agreement is obtained for Dr = 1. 

5. Comments on stellar models with Dr = 1 

Using stellar models computed with element diffusion in 
the following we give two illustrations of the use of radia- 
tive diffusivity (1) to estimate the age of Hyades and (2) 
to give properties of a sample of Pop I stellar models. 

Age and predicted abundances of the Hyades stars: So far 
we have used for the Hyades open cluster the age and 
the initial chemical composition derived from the work 
of Lebreton et al. ( |200l| ) that does not take into account 
the element diffusion. Using stellar models computed with 
Dr = 1.0, and the corrected initial helium content and 
metallicity introduced in Sect. ^, we have reestimated the 
age by fitting isochrones to the observations of Balachan- 
dran (1995). Figure || shows the isochrones. 

The inclusion of the element diffusion leads to the age 
^Hyades = 570 ± 15 Myr, i.e. an increase of w 3% with re- 
spect to the result of Lebreton et al. ( [2001 ). As already 
emphasized by Lebreton et al. ( |2001| ) one can see the sen- 
sitivity of the age to the loci of the two components of 
the 9 Tau binary system with 2 Tau located right on the 
turn-off. 
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Fig. 5. Profiles of expected surface abundances (in dex) 
with respect to effective temperature for X = He (dashed) , 
Li (full) and N (dash-dot-dot-dot-dash) for a sample of 
"Hyades" stellar models at the age of 570 My computed 
with D R = 1.0. 



Figure |^ reproduces the predictions of surface abun- 
dances for helium, lithium and nitrogen with respect to 
the effective temperature. Below T e g ^ 6 000K the exter- 
nal convective zone extends deeply enough to reach the 
area where lithium is nuclearly destroyed, resulting in a 
low surface abundance. For 6 400 ^ T e ff SS 6 600K the 
lithium abundance is close to its initial value as for the 
other elements. Above ^ 6 600K and until w 8 200K the 
depletion due to the element diffusion is more and more 
counter balanced by the mixing generated by the increas- 
ing radiative diffusivity. For stars above the turn-off, first 
the evolutionary time scale is small and the element diffu- 
sion is less efficient and second, the mixing resulting from 
the dredge up becomes more and more efficient and the 
surface abundances keep about their initial values. 

Figure |^ shows the surface abundance of lithium with 
respect to effective temperature. The so-called Li dip 
(Boesgaard & Tripicco 1986| ), located in the vicinity of 
T e g ss 6 500K is not reproduced by our models. As in 
the standard solar model, see Table [l[ the radiative dif- 
fusivity combined with element diffusion does not suffice 
to explain the lithium depletion. As clearly demonstrated 



by Charbonnel et al. (1994) the Li dip does not result 
from the gravitational settling; according to the present 



paradigm (Sills & Deliyannis 2000), it results from the 
mixing generated beneath the tachocline due to the shear 
caused by the differences of rotational velocities between 
the convection zone and the radiative interior (Talon & 
Charbonnel 1998). At the age under consideration, accord- 



ing to the rotational status of the star and to the extend 
of the external convection zone, the mixed region extends 
towards the part of the envelope where the lithium is nu- 
clearly destroyed. The other elements, except beryllium 
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Fig. 4. Surface abundances (in dex) with respect to effective temperature of C, O, Fe and Mg for a sample of "Hyades" 
stellar models at the age of 570 My computed with Dr = 0.0 (dotted), £>r = 0.5 (dashed), Z?r = 1.0 (full), Z?r = 5.0 
(dash-dot-d ash-d ot) and Z?r = 9.0 (dash-dot-dot-dot-dash). The triangles correspond to the observations of Varenne 
& Monier ( 1999 ) corrected from the metallicity of Hyades. Dashed rectangles delimit the uncertainty domains of 
observations. 



and boron, being nuclearly destroyed in the core are not 
affected by this process. The radiative diffusivity then acts 
together with the other processes to inhibit the gravita- 
tional settling, generating the mild depletion of metals as 
observed by Varenne & Monier (1999). 



For reference, Table |2j shows theoretical surface param- 
eters of our "Hyades" stellar models and reveals significant 
changes of the surface abundances according to the stellar 
mass. Therefore, at present, neither the helium mass frac- 
tion nor the metallicity have a unique value at the surface 
of the main sequence star members of the cluster. 

A sample of Pop. I stellar models: To investigate the 
change of the surface abundances resulting from the el- 
ement diffusion on the main sequence and subgiant stars 
we have computed typical evolutions of Pop. I stellar mod- 
els for respectively l.M , 1.4 M©, 1.8 M , 2.5 M Q and 



4. Mo. Models are initialized at homogeneous ZAMS with 
solar composition. Figure shows, with respect to age, the 
profiles of the surface abundance and the corresponding 
HR diagrams. Along the main sequence, the surface mass 
fractions of helium and metals are depleted, but along the 
giant branch, the surface abundances recover about their 
initial values as a consequence of the dredge-up. Therefore 
it not possible to derive a unique metallicity for all stats 
of an open cluster like Hyades. 

Note that for the model of 1 M Q at 4.65 Gyr, the age 
of the Sun, we obtain a depletion for iron of « 0.059 dex, 
consistent with the present day adopted difference be- 
tween photospheric and meteoritic observed abundances 
(Asplund et al. |2000b| ). 
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Table 2. Typical physical parameters of "Hyades" stellar models computed with Dr = 1 at age 570 Myr. M*, T e g, 
\ogg and Mboi are respectively the mass, the effective temperature (in K), the logarithm of the gravity and the 
bolometric magnitude (Mboi© = 4.75). [^], [§], [§], [§], [§] are the relative surface abundances with respect 

to the Hyades one: +0.14 dex. 
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Fig. 6. Isochrones of ZAMS (dotted), 570 Myr (full), 
585 Myr (dash-dot-dash) and 555 Myr (dashed) for the 
Hyades open cluster computed with element diffusion and 
D R = 1. 

6. Discussion 

Radiative diffusivity is a basic physical phenomenon that 
complements standard element diffusion in stellar model- 
ing. In this exploratory work, we have introduced in the 
calculation of the microscopic diffusion coefficients, an es- 
timate of the radiative diffusivity resulting from photon- 
ion collisions that acts predominantly in the external lay- 
ers of stars of high effective temperature. This simple 
treatment needs an efficiency factor, the radiative diffu- 
sivity parameter Dr. We have constrained .Dr by the fit 
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Fig. 8. Surface abundance of lithium with respect to ef- 
fective temperature of stellar models (full) of the Hyades 
computed with Dr = 1.0 at the age of 570 Myr. The la- 
bels of dots give the mass (solar unit) of models. Observed 
lithium abundances and effective temperatures of th e 
Hyades field stars (*) are taken from Balachandran ( 1995 ). 



ies presented in Sect. H strongly sugge 
to unity, Dr = 1.0 21 satisfies all th 



of observations chosen independent of rotation. The stud- 
jest that a value close 
the observational con- 
straints. We have adopted a unique value Dr. = 1 for all 
chemical element species and in all physical conditions. 

Without taking into account other physical processes 
that inhibit the gravitational sedimentation, namely mass- 
loss, or turbulent mixing generated by rotation or by some 
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Fig. 7. Panels (1) to (7) show the profiles, with respect to age, of the effective temperature, of the luminosity, of 
abundance mass fractions at the surface of a sample of Pop. I stellar models computed with = 1.0. Panel (8) shows 
the corresponding HR diagram. The labels are self explanatory. The stellar masses are respectively M+/Mq = 1.0 
(full), 1.4 (dashed), 1.8 (dot-dash-dot-dash), 2.5 (dotted), 4.0 (dash-dot-dot-dot-dash). 



other turbulent process, we find that the implementa- 
tion of the radiative diffusivity inhibits the large gravi- 
tational settling resulting from the element diffusion in 
outer layers of early-type star models with masses larger 
than ^ 1.4M . 

The radiative diffusivity is not able to solve all the 
open questions of element depletion in stellar models (e.g. 



Balachandran 1995; Sills & Deliyannis 2000 for reviews). 
Here, our goal is to emphasize on its large magnitude and 
to point out that it should be added to the other mix- 
ing processes. Radiative diffusivity simply ac ts tog ether 
with the mixing induced b y rot ation (Zahn 1992| ), the 



mass- loss (Chaboyer et al. 1999 ), the turbulent mixing 



(Scha tzman 1969|) and the radiative accelerations (Alecian 



et al. 1989) to inhibit the large element diffusion. 



The implementation of the radiative diffusivity in stel- 
lar modeling is necessary, in particular, for the calculation 
of: 



— precise models of intermediate mass stars that are tar- 
gets of asteroseismological experiments like COROT or 

EDDINGTON, 

— theoretical isochrones including element diffusion, for 
young cluster presenting stars with masses above 
1.4 M Q at the turn-off. 

Our simplified model takes into account only the mix- 
ing generated by the anisotropic component of the ra- 
diative diffusion; it ignores the radiative acceleration and 
so cannot describe the AmFm phenomenon (Michaud et 
al. h983l Alecian Il996|; Hui-Bon-Hua & Alecian (Il998l 



Richer et al. [20001 ). 

We are aware that our estimation of the radiative diffu- 
sivity is only a first approximation; owing to its efficiency 
we are confident of the fact that the diffusivity generated 
by the radiative collisions is a part of the missing physi- 
cal process responsible for the too large efficiency of the 
microscopic diffusion coefficient. Our description of the 
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radiative diffusivity is phenomenological and preliminary 
and needs further extended developments based on the 
kinetic theory of gases and radiation. 
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